Persistent homology reveals strong phylogenetic signal in 3D protein structures.

Changes that occur in proteins over time provide a phylogenetic signal that can be used to decipher their evolutionary history and the relationships between organisms. Sequence comparison is the most common way to access this phylogenetic signal, while those based on 3D structure comparisons are still in their infancy. In this study, we propose an effective approach based on Persistent Homology Theory (PH) to extract the phylogenetic information contained in protein structures. PH provides efficient and robust algorithms for extracting and comparing geometric features from noisy datasets at different spatial resolutions. PH has a growing number of applications in the life sciences, including the study of proteins (e.g. classification, folding). However, it has never been used to study the phylogenetic signal they may contain. Here, using 518 protein families, representing 22,940 protein sequences and structures, from 10 major taxonomic groups, we show that distances calculated with PH from protein structures correlate strongly with phylogenetic distances calculated from protein sequences, at both small and large evolutionary scales. We test several methods for calculating PH distances and propose some refinements to improve their relevance for addressing evolutionary questions. This work opens up new perspectives in evolutionary biology by proposing an efficient way to access the phylogenetic signal contained in protein structures, as well as future developments of topological analysis in the life sciences.

Deciphering Evolutionary Trajectories of Lactate Dehydrogenases Provides New Insights into Allostery.

Lactate dehydrogenase (LDH, EC.1.1.127) is an important enzyme engaged in the anaerobic metabolism of cells, catalyzing the conversion of pyruvate to lactate and NADH to NAD+. LDH is a relevant enzyme to investigate structure-function relationships. The present work provides the missing link in our understanding of the evolution of LDHs. This allows to explain (i) the various evolutionary origins of LDHs in eukaryotic cells and their further diversification and (ii) subtle phenotypic modifications with respect to their regulation capacity. We identified a group of cyanobacterial LDHs displaying eukaryotic-like LDH sequence features. The biochemical and structural characterization of Cyanobacterium aponinum LDH, taken as representative, unexpectedly revealed that it displays homotropic and heterotropic activation, typical of an allosteric enzyme, whereas it harbors a long N-terminal extension, a structural feature considered responsible for the lack of allosteric capacity in eukaryotic LDHs. Its crystallographic structure was solved in 2 different configurations typical of the R-active and T-inactive states encountered in allosteric LDHs. Structural comparisons coupled with our evolutionary analyses helped to identify 2 amino acid positions that could have had a major role in the attenuation and extinction of the allosteric activation in eukaryotic LDHs rather than the presence of the N-terminal extension. We tested this hypothesis by site-directed mutagenesis. The resulting C. aponinum LDH mutants displayed reduced allosteric capacity mimicking those encountered in plants and human LDHs. This study provides a new evolutionary scenario of LDHs that unifies descriptions of regulatory properties with structural and mutational patterns of these important enzymes.

Biochemical, structural and dynamical characterizations of the lactate dehydrogenase from Selenomonas ruminantium provide information about an intermediate evolutionary step prior to complete allosteric regulation acquisition in the super family of lactate and malate dehydrogenases.

In this work, we investigated the lactate dehydrogenase (LDH) from Selenomonas ruminantium (S. rum), an enzyme that differs at key amino acid positions from canonical allosteric LDHs. The wild type (Wt) of this enzyme recognises pyuvate as all LDHs. However, introducing a single point mutation in the active site loop (I85R) allows S. Rum LDH to recognize the oxaloacetate substrate as a typical malate dehydrogenase (MalDH), whilst maintaining homotropic activation as an LDH. We report the tertiary structure of the Wt and I85RLDH mutant. The Wt S. rum enzyme structure binds NADH and malonate, whilst also resembling the typical compact R-active state of canonical LDHs. The structure of the mutant with I85R was solved in the Apo State (without ligand), and shows no large conformational reorganization such as that observed with canonical allosteric LDHs in Apo state. This is due to a local structural feature typical of S. rum LDH that prevents large-scale conformational reorganization. The S. rum LDH was also studied using Molecular Dynamics simulations, probing specific local deformations of the active site that allow the S. rum LDH to sample the T-inactive state. We propose that, with respect to the LDH/MalDH superfamily, the S. rum enzyme possesses a specificstructural and dynamical way to ensure homotropic activation.

Protein Conformational Space at the Edge of Allostery: Turning a Nonallosteric Malate Dehydrogenase into an "Allosterized" Enzyme Using Evolution-Guided Punctual Mutations.

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

Phylogenetics and biochemistry elucidate the evolutionary link between l-malate and l-lactate dehydrogenases and disclose an intermediate group of sequences with mix functional properties.

The NAD(P)-dependent malate dehydrogenases (MDH) (EC 1.1.1.37) and NAD-dependent lactate dehydrogenases (LDH) (EC. 1.1.1.27) form a large superfamily that has been characterized in organisms belonging to the three Domains of Life. MDH catalyzes the reversible conversion of the oxaloacetate into malate, while LDH operates at the late stage of glycolysis by converting pyruvate into lactate. Phylogenetic studies proposed that the LDH/MDH superfamily encompasses five main groups of enzymes. Here, starting from 16,052 reference proteomes, we reinvestigated the relationships between MDH and LDH. We showed that the LDH/MDH superfamily encompasses three main families: MDH1, MDH2, and a large family encompassing MDH3, LDH, and L-2-hydroxyisocaproate dehydrogenases (HicDH) sequences. An in-depth analysis of the phylogeny of the MDH3/LDH/HicDH family and of the nature of three important amino acids, located within the catalytic site and involved in binding and substrate discrimination, revealed a large group of sequences displaying unexpected combinations of amino acids at these three critical positions. This group branched in-between canonical MDH3 and LDH sequences. The functional characterization of several enzymes from this intermediate group disclosed a mix of functional properties, indicating that the MDH3/LDH/HicDH family is much more diverse than previously thought, and blurred the frontier between MDH3 and LDH enzymes. Present-days enzymes of the intermediate group are a valuable material to study the evolutionary steps that led to functional diversity and emergence of allosteric regulation within the LDH/MDH superfamily.

Biochemical, structural and dynamical studies reveal strong differences in the thermal-dependent allosteric behavior of two extremophilic lactate dehydrogenases.

In this work, we combined biochemical and structural investigations with molecular dynamics (MD) simulations to analyze the very different thermal-dependent allosteric behavior of two lactate dehydrogenases (LDH) from thermophilic bacteria. We found that the enzyme from Petrotoga mobilis (P. mob) necessitates an absolute requirement of the allosteric effector (fructose 1, 6-bisphosphate) to ensure functionality. In contrast, even without allosteric effector, the LDH from Thermus thermophilus (T. the) is functional when the temperature is raised. We report the crystal structure of P. mob LDH in the Apo state solved at 1.9 Å resolution. We used this structure and the one from T. the, obtained previously, as a starting point for MD simulations at various temperatures. We found clear differences between the thermal dynamics, which accounts for the behavior of the two enzymes. Our work demonstrates that, within an allosteric enzyme, some areas act as local gatekeepers of signal transmission, allowing the enzyme to populate either the T-inactive or the R-active states with different degrees of stringency.

Resurrection of Ancestral Malate Dehydrogenases Reveals the Evolutionary History of Halobacterial Proteins: Deciphering Gene Trajectories and Changes in Biochemical Properties.

Extreme halophilic Archaea thrive in high salt, where, through proteomic adaptation, they cope with the strong osmolarity and extreme ionic conditions of their environment. In spite of wide fundamental interest, however, studies providing insights into this adaptation are scarce, because of practical difficulties inherent to the purification and characterization of halophilic enzymes. In this work, we describe the evolutionary history of malate dehydrogenases (MalDH) within Halobacteria (a class of the Euryarchaeota phylum). We resurrected nine ancestors along the inferred halobacterial MalDH phylogeny, including the Last Common Ancestral MalDH of Halobacteria (LCAHa) and compared their biochemical properties with those of five modern halobacterial MalDHs. We monitored the stability of these various MalDHs, their oligomeric states and enzymatic properties, as a function of concentration for different salts in the solvent. We found that a variety of evolutionary processes, such as amino acid replacement, gene duplication, loss of MalDH gene and replacement owing to horizontal transfer resulted in significant differences in solubility, stability and catalytic properties between these enzymes in the three Halobacteriales, Haloferacales, and Natrialbales orders since the LCAHa MalDH. We also showed how a stability trade-off might favor the emergence of new properties during adaptation to diverse environmental conditions. Altogether, our results suggest a new view of halophilic protein adaptation in Archaea.

Protein thermal stability.

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Remote oxidative modifications induced by oxygen free radicals modify T/R allosteric equilibrium of a hyperthermophilic lactate dehydrogenase.

L-Lactate dehydrogenase (LDH) is a model protein allowing to shed light on the fundamental molecular mechanisms that drive the acquisition, evolution and regulation of enzyme properties. In this study, we test the hypothesis of a link between thermal stability of LDHs and their capacity against unfolding induced by reactive oxygen species (ROS) generated by γ-rays irradiation. By using circular dichroism spectroscopy, we analysed that high thermal stability of a thermophilic LDH favours strong resistance against ROS-induced unfolding, in contrast to its psychrophilic and mesophilic counterparts that are less resistant. We suggest that a protein's phenotype linking strong thermal stability and resistance against ROS damages would have been a selective evolutionary advantage. We also find that the enzymatic activity of the thermophilic LDH that is strongly resistant against ROS-unfolding is very sensitive to inactivation by irradiation. To address this counter-intuitive observation, we combined mass spectrometry analyses and enzymatic activity measurements. We demonstrate that the dramatic change on LDH activity was linked to remote chemical modifications away from the active site, that change the equilibrium between low-affinity tense (T-inactive) and high-affinity relaxed (R-active) forms. We found the T-inactive thermophilic enzyme obtained after irradiation can recover its LDH activity by addition of the allosteric effector 1, 6 fructose bis phosphate. We analyse our data within the general framework of allosteric regulation, which requires that an enzyme in solution populates a large diversity of dynamically-interchanging conformations. Our work demonstrates that the radiation-induced inactivation of an enzyme is controlled by its dynamical properties.

The archaeal LDH-like malate dehydrogenase from Ignicoccus islandicus displays dual substrate recognition, hidden allostery and a non-canonical tetrameric oligomeric organization.

The NAD(P)-dependent malate dehydrogenases (MalDHs) and NAD-dependent lactate dehydrogenases (LDHs) are homologous enzymes involved in central metabolism. They display a common protein fold and the same catalytic mechanism, yet have a stringent capacity to discriminate between their respective substrates. The MalDH/LDH superfamily is divided into several phylogenetically related groups. It has been shown that the canonical LDHs and LDH-like group of MalDHs are primarily tetrameric enzymes that diverged from a common ancestor. In order to gain understanding of the evolutionary history of the LDHs and MalDHs, the biochemical properties and crystallographic structure of the LDH-like MalDH from the hyperthermophilic archaeon Ignicoccus islandicus (I. isl) were determined. I. isl MalDH recognizes oxaloacetate as main substrate, but it is also able to use pyruvate. Surprisingly, with pyruvate, the enzymatic activity profile looks like that of allosteric LDHs, suggesting a hidden allosteric capacity in a MalDH. The I. isl MalDH tetrameric structure in the apo state is considerably different from those of canonical LDH-like MalDHs and LDHs, representing an alternative oligomeric organization. A comparison with MalDH and LDH counterparts provides strong evidence that the divergence between allosteric and non-allosteric members of the superfamily involves homologs with intermediate, atypical properties.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase.

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Crystallophore: a versatile lanthanide complex for protein crystallography combining nucleating effects, phasing properties, and luminescence.

Macromolecular crystallography suffers from two major issues: getting well-diffracting crystals and solving the phase problem inherent to large macromolecules. Here, we describe the first example of a lanthanide complex family named "crystallophore" (Xo4), which contributes to tackling both bottlenecks. This terbium complex, Tb-Xo4, is an appealing agent for biocrystallography, combining the exceptional phasing power of the Tb(iii) heavy atom with powerful nucleating properties, providing ready-to-use crystals for structure determination. Furthermore, protein/Tb-Xo4 co-crystals can be easily detected and discriminated from other crystalline by-products using luminescence. We demonstrate the potential of this additive for the crystallisation and structure determination of eight proteins, two of whose structures were unknown.

SAXS/SANS on Supercharged Proteins Reveals Residue-Specific Modifications of the Hydration Shell.

Water molecules in the immediate vicinity of biomacromolecules, including proteins, constitute a hydration layer characterized by physicochemical properties different from those of bulk water and play a vital role in the activity and stability of these structures, as well as in intermolecular interactions. Previous studies using solution scattering, crystallography, and molecular dynamics simulations have provided valuable information about the properties of these hydration shells, including modifications in density and ionic concentration. Small-angle scattering of x-rays (SAXS) and neutrons (SANS) are particularly useful and complementary techniques to study biomacromolecular hydration shells due to their sensitivity to electronic and nuclear scattering-length density fluctuations, respectively. Although several sophisticated SAXS/SANS programs have been developed recently, the impact of physicochemical surface properties on the hydration layer remains controversial, and systematic experimental data from individual biomacromolecular systems are scarce. Here, we address the impact of physicochemical surface properties on the hydration shell by a systematic SAXS/SANS study using three mutants of a single protein, green fluorescent protein (GFP), with highly variable net charge (+36, -6, and -29). The combined analysis of our data shows that the hydration shell is locally denser in the vicinity of acidic surface residues, whereas basic and hydrophilic/hydrophobic residues only mildly modify its density. Moreover, the data demonstrate that the density modifications result from the combined effect of residue-specific recruitment of ions from the bulk in combination with water structural rearrangements in their vicinity. Finally, we find that the specific surface-charge distributions of the different GFP mutants modulate the conformational space of flexible parts of the protein.

New insights into the mechanism of substrates trafficking in Glyoxylate/Hydroxypyruvate reductases.

Glyoxylate accumulation within cells is highly toxic. In humans, it is associated with hyperoxaluria type 2 (PH2) leading to renal failure. The glyoxylate content within cells is regulated by the NADPH/NADH dependent glyoxylate/hydroxypyruvate reductases (GRHPR). These are highly conserved enzymes with a dual activity as they are able to reduce glyoxylate to glycolate and to convert hydroxypyruvate into D-glycerate. Despite the determination of high-resolution X-ray structures, the substrate recognition mode of this class of enzymes remains unclear. We determined the structure at 2.0 Å resolution of a thermostable GRHPR from Archaea as a ternary complex in the presence of D-glycerate and NADPH. This shows a binding mode conserved between human and archeal enzymes. We also determined the first structure of GRHPR in presence of glyoxylate at 1.40 Å resolution. This revealed the pivotal role of Leu53 and Trp138 in substrate trafficking. These residues act as gatekeepers at the entrance of a tunnel connecting the active site to protein surface. Taken together, these results allowed us to propose a general model for GRHPR mode of action.

Xanthan Exopolysaccharide: Cu(2+) Complexes Affected from the pH-Dependent Conformational State; Implications for Environmentally Relevant Biopolymers.

The conformational impact of environmental biopolymers on metal sorption was studied through Cu sorption on xanthan. The apparent Cu(2+) complexation constant (logK; Cu(2+) + L(-) ↔ CuL(+)) decreased from 2.9 ± 0.1 at pH 3.5 to 2.5 ± 0.1 at pH 5.5 (ionic strength I = 0.1). This behavior is in apparent contradiction with basic thermodynamics, as usually the higher the pH the more cations bind. Our combined titration, circular dichroism and dynamic light scattering study indicated that the change observed in Cu bond strength relates to a conformational change of the structure of xanthan, which generates more chelating sites at pH 3.5 than at pH 5.5. This hypothesis was validated by the fact that the Cu sorption constants on xanthan were always higher than those measured on a mixture of pyruvic and glucuronic acids (logK = 2.2), which are the two constitutive ligands present in the xanthan monomer. This study shows the role of the structural conformation of natural biopolymers in metal bond strength. This finding may help to better predict the fate of Cu and other metals in acidic environmental settings such as aquatic media affected by acid mine drainage, as well as peats and acidic soils, and to better define optimal conditions for bioremediation processes.

Interface matters: the stiffness route to stability of a thermophilic tetrameric malate dehydrogenase.

In this work we investigate by computational means the behavior of two orthologous bacterial proteins, a mesophilic and a thermophilic tetrameric malate dehydrogenase (MalDH), at different temperatures. Namely, we quantify how protein mechanical rigidity at different length- and time-scales correlates to protein thermophilicity as commonly believed. In particular by using a clustering analysis strategy to explore the conformational space of the folded proteins, we show that at ambient conditions and at the molecular length-scale the thermophilic variant is indeed more rigid that the mesophilic one. This rigidification is the result of more efficient inter-domain interactions, the strength of which is further quantified via ad hoc free energy calculations. When considered isolated, the thermophilic domain is indeed more flexible than the respective mesophilic one. Upon oligomerization, the induced stiffening of the thermophilic protein propagates from the interface to the active site where the loop, controlling the access to the catalytic pocket, anchors down via an extended network of ion-pairs. On the contrary in the mesophilic tetramer the loop is highly mobile. Simulations at high temperature, could not re-activate the mobility of the loop in the thermophile. This finding opens questions on the similarities of the binding processes for these two homologues at their optimal working temperature and suggests for the thermophilic variant a possible cooperative role of cofactor/substrate.

Phylogenetically driven sequencing of extremely halophilic archaea reveals strategies for static and dynamic osmo-response.

Organisms across the tree of life use a variety of mechanisms to respond to stress-inducing fluctuations in osmotic conditions. Cellular response mechanisms and phenotypes associated with osmoadaptation also play important roles in bacterial virulence, human health, agricultural production and many other biological systems. To improve understanding of osmoadaptive strategies, we have generated 59 high-quality draft genomes for the haloarchaea (a euryarchaeal clade whose members thrive in hypersaline environments and routinely experience drastic changes in environmental salinity) and analyzed these new genomes in combination with those from 21 previously sequenced haloarchaeal isolates. We propose a generalized model for haloarchaeal management of cytoplasmic osmolarity in response to osmotic shifts, where potassium accumulation and sodium expulsion during osmotic upshock are accomplished via secondary transport using the proton gradient as an energy source, and potassium loss during downshock is via a combination of secondary transport and non-specific ion loss through mechanosensitive channels. We also propose new mechanisms for magnesium and chloride accumulation. We describe the expansion and differentiation of haloarchaeal general transcription factor families, including two novel expansions of the TATA-binding protein family, and discuss their potential for enabling rapid adaptation to environmental fluxes. We challenge a recent high-profile proposal regarding the evolutionary origins of the haloarchaea by showing that inclusion of additional genomes significantly reduces support for a proposed large-scale horizontal gene transfer into the ancestral haloarchaeon from the bacterial domain. The combination of broad (17 genera) and deep (≥5 species in four genera) sampling of a phenotypically unified clade has enabled us to uncover both highly conserved and specialized features of osmoadaptation. Finally, we demonstrate the broad utility of such datasets, for metagenomics, improvements to automated gene annotation and investigations of evolutionary processes.

An experimental point of view on hydration/solvation in halophilic proteins.

Protein-solvent interactions govern the behaviors of proteins isolated from extreme halophiles. In this work, we compared the solvent envelopes of two orthologous tetrameric malate dehydrogenases (MalDHs) from halophilic and non-halophilic bacteria. The crystal structure of the MalDH from the non-halophilic bacterium Chloroflexus aurantiacus (Ca MalDH) solved, de novo, at 1.7 Å resolution exhibits numerous water molecules in its solvation shell. We observed that a large number of these water molecules are arranged in pentagonal polygons in the first hydration shell of Ca MalDH. Some of them are clustered in large networks, which cover non-polar amino acid surface. The crystal structure of MalDH from the extreme halophilic bacterium Salinibacter ruber (Sr) solved at 1.55 Å resolution shows that its surface is strongly enriched in acidic amino acids. The structural comparison of these two models is the first direct observation of the relative impact of acidic surface enrichment on the water structure organization between a halophilic protein and its non-adapted counterpart. The data show that surface acidic amino acids disrupt pentagonal water networks in the hydration shell. These crystallographic observations are discussed with respect to halophilic protein behaviors in solution.

Sampling the conformational energy landscape of a hyperthermophilic protein by engineering key substitutions.

Proteins exist as a dynamic ensemble of interconverting substates, which defines their conformational energy landscapes. Recent work has indicated that mutations that shift the balance between conformational substates (CSs) are one of the main mechanisms by which proteins evolve new functions. In the present study, we probe this assertion by examining phenotypic protein adaptation to extreme conditions, using the allosteric tetrameric lactate dehydrogenase (LDH) from the hyperthermophilic bacterium Thermus thermophilus (Tt) as a model enzyme. In the presence of fructose 1, 6 bis-phosphate (FBP), allosteric LDHs catalyze the conversion of pyruvate to lactate with concomitant oxidation of nicotinamide adenine dinucleotide, reduced form (NADH). The catalysis involves a structural transition between a low-affinity inactive "T-state" and a high-affinity active "R-state" with bound FBP. During this structural transition, two important residues undergo changes in their side chain conformations. These are R171 and H188, which are involved in substrate and FBP binding, respectively. We designed two mutants of Tt-LDH with one ("1-Mut") and five ("5-Mut") mutations distant from the active site and characterized their catalytic, dynamical, and structural properties. In 1-Mut Tt-LDH, without FBP, the K(m)(Pyr) is reduced compared with that of the wild type, which is consistent with a complete shifting of the CS equilibrium of H188 to that observed in the R-state. By contrast, the CS populations of R171, k(cat) and protein stability are little changed. In 5-Mut Tt-LDH, without FBP, K(m)(Pyr) approaches the values it has with FBP and becomes almost temperature independent, k(cat) increases substantially, and the CS populations of R171 shift toward those of the R-state. These changes are accompanied by a decrease in protein stability at higher temperature, which is consistent with an increased flexibility at lower temperature. Together, these results show that the thermal properties of an enzyme can be strongly modified by only a few or even a single mutation, which serve to alter the equilibrium and, hence, the relative populations of functionally important native-state CSs, without changing the nature of the CSs themselves. They also provide insights into the types of mutational pathways by which protein adaptation to temperature is achieved.